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ABSTRACT 

We estimate the acceleration on the Local Group (LG) from the Two Micron All Sky Redshift 
Survey (2MRS). The sample used includes about 23,200 galaxies with extinction corrected 
magnitudes brighter than K s = 11.25 and it allows us to calculate the flux weighted dipole. 
The near-infrared flux weighted dipoles are very robust because they closely approximate a 
mass weighted dipole, bypassing the effects of redshift distortions and require no preferred 
reference frame. This is combined with the redshift information to determine the change in 
dipole with distance. The misalignment angle between the LG and the CMB dipole drops 
to 12°±7° at around 50 ft, Mpc, but then increases at larger distances, reaching 21°±8° 
at around 130 ft. Mpc. Exclusion of the galaxies Maffei 1, Maffei 2, Dwingeloo 1, IC342 
and M87 brings the resultant flux dipole to 14°±7° away from the CMB velocity dipole 
In both cases, the dipole seemingly converges by 60 /i _1 Mpc. Assuming convergence, the 
comparison of the 2MRS flux dipole and the CMB dipole provides a value for the combination 
of the mass density and luminosity bias parameters £! m 0,6 /£>L = 0.40 ± 0.09. 

Key words: methods:data analysis- cosmology: observations - large-scale structure of uni- 
verse - galaxies: Local Group - infrared:galaxies 



1 INTRODUCTION 

The most popular mechanism for the formation of large-scale struc- 
ture and motions in the Universe is the gravitational growth of pri- 
mordial density perturbations. According to this paradigm, if the 
density perturbations are small enough to be approximated by a 
linear theory, then the peculiar acceleration vector g(r) induced by 
the matter distribution around position r is related to the mass by 



oo 

?(r)=Gp J dVMOp: 



(1) 



where p is the mean matter density and 5 m (r) = (p m (r) — p)/p is 
the density contrast of the mass perturbations. In linear theory, the 



peculiar velocity field, v(r), is proportional to the peculiar acceler- 
ation: 

where Ho = 100 h kms~ 1 Mpc~ 1 is the Hubble constant and 
/(fi m ) ~ ^m 6 is the logarithmic derivative of the amplitude of 
the growing mode of the perturbations in mass with respect to the 
scale factor (Peebles 1980). The factor /(fi m ) is only weakly de- 
pendent on the cosmological constant (Lahav et al, 1991). 

During the past twenty five years, particular attention has been 
paid to the study of the gravitational acceleration and the peculiar 
velocity vectors of the Local Group (LG) of galaxies. It is now 
widely accepted that the cosmic microwave background (CMB) 
dipole is a Doppler effect arising from the motion of the Sun (but 
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see e.g. Gunn 1988 and Paczyriski & Piran 1990 who argue that 
the CMB dipole is of primordial origin). In this case, the dipole 
anisotropy of the CMB is a direct and accurate measurement of 
the LG peculiar velocity (c.f. Conklin 1969 and Henry 1971). The 
LG acceleration can also be estimated using surveys of the galax- 
ies tracing the density inhomogeneities responsible for the accel- 
eration. By comparing the CMB velocity vector with the acceler- 
ation vector 1 obtained from the galaxy surveys, it is possible to 
investigate the cause of the LG motion and its cosmological im- 
plications. This technique was first applied by Yahil, Sandage & 
Tammann (1980) using the Revised Shapley-Ames catalogue and 
later by Davis & Huchra (1982) using the CfA catalogue. Both cat- 
alogues were two-dimensional and the analyses were done using 
galaxy fluxes. Since both the gravity and the flux are inversely pro- 
portional to the square of the distance, the dipole vector can be cal- 
culated by summing the flux vectors and assuming an average value 
for the mass-to-light ratio. Lahav (1987) applied the same method 
to calculate the dipole anisotropy using maps based on three galaxy 
catalogues, UGC, ESO, MCG. The most recent application of the 
galaxy flux dipole analysis was carried out by Mailer et al. (2003), 
using the two-dimensional Two Micron All-Sky Survey (2MASS) 
extended source catalogue (XSC), with a limiting magnitude of 
K s = 13.57. They found that the LG dipole direction is 16° away 
that of the CMB. 

Our ability to study the LG motion was greatly advanced by 
the whole-sky galaxy samples derived from IRAS Galaxy Cata- 
logues. Yahil, Walker & Rowan-Robinson (1986), Meiksin & Davis 

(1986) , Harmon, Lahav & Meurs (1987), Villumsen & Strauss 

(1987) and Lahav, Rowan-Robinson & Lynden-Bell (1988) used 
the two-dimensional IRAS data to obtain the LG dipole. The 
dipole vectors derived by these authors are in agreement with each 
other and the CMB dipole vector to within 10°-30° degrees. The 
inclusion of galaxy redshifts in the dipole analyses allowed the 
estimation of the distance at which most of the peculiar veloc- 
ity of the LG is generated (the convergence depth). However, the 
estimates of the convergence depth from various data sets have 
not agreed. Strauss et al. (1992, IRAS sample), Webster, Lahav 
& Fisher (1997, IRAS sample), Lynden-Bell, Lahav & Burstein 
(1989, optical sample) and da Costa et al. (2000, a sample of early- 
type galaxies) suggested that the LG acceleration is mostly due to 
galaxies < 50/i~ 1 Mpc, while other authors such as Scaramella, 
Vettolani & Zamorani (1994, Abell/ACO cluster sample), Bran- 
chini & Plionis (1996, Abell/ACO cluster sample), Kocevski et al. 
(2004) and Kocevski & Ebeling (2005, both using samples of X-ray 
clusters) claimed that there is significant contribution to the dipole 
from depths of up to ~200 ft _1 Mpc. 

Dipole analyses are often used to estimate the combination of 
matter density and biasing parameters f2 m and b. In theory, one can 
equate the velocity inferred from the CMB measurements with the 
value derived from a galaxy survey and obtain a value for /3. In 
practice, however, the galaxy surveys do not measure the true total 
velocity due their finite depth (e.g. Lahav, Kaiser & Hoffman 1990 
and Juszkiewicz, Vittorio & Wyse 1990). The true vlg, as obtained 
from the CMB dipole, arises from structure on all scales including 
structures further away than the distance a galaxy survey can ac- 
curately measure. Furthermore, the magnitude/flux/diameter limit 
of the survey and any completeness variations over the sky intro- 



1 Both the CMB velocity and the gravitational acceleration on the LG have 
units of velocity and are commonly referred to as 'dipoles'. Hereafter, the 
terms 'LG velocity' and 'LG dipole' will be used interchangeably. 



duce selection effects and biases to the calculations. These effects 
amplify the errors at large distances (for redshift surveys) and faint 
magnitudes (for two dimensional surveys) where the sampling of 
the galaxy distribution becomes more sparse. This, combined with 
the fact that we sample discretely from an underlying continuous 
mass distribution leads to an increase in shot noise error. There 
may also be a significant contribution to the dipole from galaxies 
behind the Galactic Plane (the zone of avoidance). The analysis of 
the convergence of the dipole is further complicated by the redshift 
distortions on small and large scales which introduce systematic er- 
rors to the derived dipole (the rocket effect, Kaiser 1987). The fol- 
lowing sections discuss these effects in the context of two different 
models of biasing. 

In this paper, we use the Two Micron All-Sky Redshift Survey 
(2MRS, Huchra etal. 2005) 2 to study the LG dipole. The inclusion 
of the redshift data allows the calculation of the selection effects of 
the survey as a function of distance and enables the study of conver- 
gence and thus improves the analysis of Mailer et al. (2003). The 
paper is structured as follows: The Two Micron Redshift Survey is 
described in Section|2| Section 3 discusses the method used in the 
analysis including the different weighting schemes, the rocket ef- 
fect and the choice of reference frames. The results are presented 
in Section The final section includes some concluding remarks 
and plans for future work. 



2 THE TWO MICRON ALL-SKY REDSHIFT SURVEY 

The Two Micron All-Sky Redshift Survey (2MRS) is the densest 
all-sky redshift survey to date. The galaxies in the northern ce- 
lestial hemisphere are being observed mainly by the FLWO 1.5-m 
telescope and at low latitudes by the CTIO. In the southern hemi- 
sphere, most galaxies are observed as a part of the six degree field 
galaxy survey (6dFGS, Jones et al. 2004) conducted by the An- 
glo Australian Observatory. The first phase of the 2MRS is now 
completed. In this phase we obtained redshifts for approximately 
23,000 2MASS galaxies from a total sample of about 24,800 galax- 
ies with extinction corrected magnitudes (Schlegel, Finkbeiner & 
Davis 1998) brighter than K s = 11.25. This magnitude limit cor- 
responds to a median redshift of z ~ 0.02 (~ 60ft" 1 Mpc). The 
majority of the 1600 galaxies that remain without redshifts are at 
very low galactic latitudes or obscured/confused by the dust and 
the high stellar density towards the Galactic Centre. FigureQshows 
all the objects in the 2MRS in Galactic Aitoff Projection. Galaxies 
with z < 0.01 are plotted in red, 0.01 < z ^ 0.025 are plotted in 
blue, 0.025 < z < 0.05 are plotted in green and z ^ 0.05 are plot- 
ted in magenta. Galaxies without measured redshifts are plotted in 
black. The 2MRS can be compared with the deeper 2MASS galaxy 
catalogue (K<14th mag) shown in Jarrett (2004, Figure 1). 

2.1 Survey Completeness 

The 2MASS 3 has great photometric uniformity and an unprece- 
dented integral sky coverage. The photometric uniformity is better 
than 4% over the sky including the celestial poles (e.g. Jarrett et al. 
2000a, 2003). The uniform completeness of the galaxy sample is 

2 This work is based on observations made at the Cerro Tololo Interameri- 
can Observatory (CTIO), operated for the US National Science Foundation 
by the Association of Universities for Research in Astronomy. 

3 The 2MASS database and the full documentation are available on the 
WWW at http//www.ipac. caltech.edu/2mass. 
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Figure 1. All Objects in the 2MASS Redshift Catalogue in Galactic Aitoff Projection. Galaxies with z ^ 0.01 are plotted in red, 0.01 < z ^ 0.025 are 
plotted in blue, 0.025 < z < 0.05 are plotted in green and z 0.05 are plotted in magenta. Galaxies without measured redshifts are plotted in black. The 
masked region is outlined by dashed lines. 



limited by the presence of the foreground stars. For a typical high 
latitude sky less than 2% of the area is masked by stars. These miss- 
ing regions are accounted for using a coverage map, defined as the 
fraction of the area of an 8' x 8'pixel that is not obscured by stars 
brighter than 10th mag. Galaxies are then weighted by the inverse 
of the completeness although the analysis is almost unaffected by 
this process as the completeness ratio is very close to one for most 
parts of the sky. 

The stellar contamination of the catalogue is low and is re- 
duced further by manually inspecting the objects below a redshift 
of cz — 200 km s -1 . The foreground stellar confusion is highest 
at low Galactic latitudes, resulting in decreasing overall complete- 
ness of the 2MASS catalogue (e.g. Jarrett et al. 2000b) and conse- 
quently the 2MRS sample 4 . Stellar confusion also produces colour 
bias in the 2MASS galaxy photometry (Cambresy, Jarrett & Be- 
ichman 2005) but this bias should not be significant for the 2MRS 
because of its relatively bright magnitude limit. 

In order to account for incompleteness at low Galactic lati- 
tudes we fill the Zone of Avoidance (the plane where |6| < 5° and 
|6| < 10° in the region |Z| < 30°) with galaxies. We keep the 



4 See Mailer et al. (2005) who reduce the stellar contamination in the 
2MASS XSC by cross-correlating stars with galaxy density. 



galaxies with observed redshifts and apply two different methods 
to compensate for the unobserved (masked) sky: 

• Method 1: The masked region is filled with galaxies whose 
fluxes and redshifts are chosen randomly from the whole data set. 
These galaxies are placed at random locations within the masked 
area. The masked region has the same average density of galaxies 
as the rest of the sky. 

• Method 2: The masked region is filled following Yahil et al. 
(1991). The area is divided into 36 bins of 10° in longitude. In 
each angular bin, the distance is divided into bins of 1000 kms -1 . 
The galaxies in each longitude/distance bin are then sampled from 
the corresponding longitude/distance bins in the adjacent strips 
-\bmasked\ - 10° < 6 < \b maskcd \ + 10° (where \b maaked \ = 5° 
or |& m asfced| = 10°). These galaxies are then placed in random lat- 
itudes within the mask region. This procedure gives similar results 
to the more elaborate method of Wiener reconstruction across the 
zone of avoidance (Lahav et al. 1994). The number of galaxies in 
each masked bin is set to a random Poisson deviate whose mean 
equals to the mean number of galaxies in the adjacent unmasked 
strips. This procedure is carried out to mimic the shot noise effects. 

In reality, the shape of the Zone of Avoidance is not as symmetric as 
defined in this paper with Galactic Bulge centred at I w +5° and 
with latitude offsets (see Kraan-Korteweg 2005). However, since 
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Figure 2. Top: Objects in the 2MASS Redshift Catalogue used in the analysis in a Galactic Aitoff Projection, including the random galaxies generated by 
using the first technique. Galaxies with z 0.01 are plotted in red, 0.01 < z sC 0.025 are plotted in blue, 0.025 < z < 0.05 are plotted in green and 
z J? 0.05 are plotted in magenta. Bottom: Same as the top plot but now including the random galaxies generated in using the second technique. The regions 
inside the dashed lines in the plots are masked out and replaced with the random galaxies shown in the plots. There are 21510 galaxies in each plot. 



we keep the galaxies in the masked regions, our dipole determi- 
nations should not be greatly influenced by assuming a symmetric 
mask. We test this by changing the centre of the Galactic bulge. We 
confirm that our results are not affected. Figure[2|shows the 2MRS 
galaxies used in the analyses in a Galactic Aitoff projection. The 
galaxies in masked regions are generated using the first (top plot) 
and the second method (bottom plot). 



2.2 Magnitude and Flux Conversions 

The 2MRS uses the 2MASS magnitude K 20 , which is defined 5 
as the magnitude inside the circular isophote corresponding to a 
surface brightness of (ak s = 20mag arcsec - (e.g. Jarrett et al. 



2000a). The isophotal magnitudes underestimate the total luminos- 
ity by 10% for the early-type and 20% for the late-type galaxies 
(Jarrett et al. 2003). Following Kochanek et al. (2001, Appendix), 
the offset of A = —0.20 ± 0.04 is added to the K 2 o magnitudes. 
The galaxy magnitudes are corrected for Galactic extinction us- 
ing the dust maps of Schlegel, Finkbeiner & Davis (1998) and an 
extinction correction coefficient of Rk = 0.35 (Cardelli, Clay- 
ton & Mathis 1989). As expected, the extinction corrections are 
small for the 2MRS sample. The K s band fc-correction is derived 
by Kochanek et al. (2001) based on the stellar population models 
of Worthey (1994). The k-correction of k(z) = -6.01og(l + z), 
is independent of galaxy type and valid for z < 0.25. 

The fluxes S are computed from the apparent magnitudes us- 
ing 



Column 17 (k_m_k20fc) in the 2MASS XSC S = 5(0 mag) W -" a ( k 2b+zpO) (3) 
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2.5 Assigning distances to nearby Galaxies 

In order to reduce the distance conversion errors, we cross- 
identified 35 galaxies which have HST Key Project distances (see 
Freedman et al. 2001 and the references therein) and 110 galaxies 
which have distance measurements compiled from several sources 
(see Karachentsev et al. 2004 and the references therein). We as- 
sign these galaxies measured distances instead of converting them 
from redshifts. In addition, we identify nine blue-shifted galaxies in 
the 2MRS which are members of the Virgo cluster. We assign these 
galaxies the distance to the centre of Virgo (~ 15.4 Mpc, Freedman 
et al. 2001). Finally, there are four remaining blue-shifted galaxies 
without known distance measurements which are assigned to 1.18 
h^ 1 Mpc 6 . Thus by assigning distances to galaxies, we do not 
need to exclude any non-LG galaxy from the analysis. 



Figure 3. Redshift histogram for 2MRS galaxies and a least squares fit 
(Equation[4j to the data (black). For comparison, also plotted is a redshift 
histogram for PSCz galaxies (Saunders et al. 2000) (red). 



where the zero point offset is ZPO = 0.017 ± 0.005 and 
S^Omag) = 1.122 x 10" 14 ± 1.891 x 10" 16 Wcm- 2 for the 
K B band (Cohen, Wheaton & Megeath 2003). 



2.3 The Redshift Distribution and the Selection Function 

The redshift distribution of the 2MRS is shown in Figure [3] The 
IRAS PSCz survey redshift distribution (Saunders et al. 2000) is 
also plotted for comparison. The 2MRS samples the galaxy distri- 
bution better than the PSCz survey out to cz = 15000 kms -1 . The 
selection function of the survey (i.e. the probability of detecting a 
galaxy as a function of distance) is modeled using a parametrised 
fit to the redshift distribution: 



dN(z) = Az 1 exp 



dz 



(4) 



with best-fit parameters of A = 116000 ± 4000, a = 2.108 ± 
0.003, 7 = 1.125±0.025 and z c = 0.025±0.001. This best-fit is 
also shown in Figure [3] (solid line). The overall selection function 
(j)(r) is the redshift distribution divided by the volume element 



n s r 2 



dN 

dz 



(5) 



where f2 s (w 4-7T steradians) is the solid angle of the survey and r is 
the comoving distance. 



2.4 Taking out the Local Group Galaxies 

Galaxies that are members of the Local Group need to be removed 
from the 2MRS catalogue to maintain the internal consistency of 
the analysis. We used the Local Group member list of thirty five 
galaxies (including Milky Way) given in Courteau and Van den 
Bergh (1999) to identify and remove eight LG members (IC 10, 
NGC 147, NGC 185, NGC 205, NGC 6822, M31, M32, M33) from 
our analysis. In Section 4, we will calculate the acceleration on to 
the Milky Way due to these LG members. 



3 THE METHODS AND WEIGHTING SCHEMES 

In order to compare the CMB and the LG dipoles, it is necessary 
to postulate a relation between the galaxy distribution and the un- 
derlying mass distribution. In this paper, we will use both a number 
weighted and a. flux weighted prescription. Although quite similar 
in formulation, these schemes are based on very different models 
for galaxy formation. The number weighted prescription assumes 
that the mass distribution in the Universe is a continuous density 
field and that the galaxies sample this field in a Poisson way. On 
the other hand, the flux weighted model is based on the assumption 
that the mass in the Universe is entirely locked to the mass of the 
halos of the luminous galaxies. 



3.1 Number Weighted Dipole 

It is commonly assumed that the galaxy and the mass distributions 
in the Universe are directly proportional to each other and are re- 
lated by a proportionality constant 7 , the linear bias parameter b: 
8n/n = S B — bS m . In this case, Equation [2] for the LG can be 
rewritten as 



VLG 



H p 
4ir 



d 3 r'5 



(6) 



where p = tt m °- 6 /b. 

For the number weighted model, incomplete sampling due to 
the magnitude limit is described by the selection function, <j){r) 
given in Section l2~3l Each galaxy i is assigned a weight: 



(7) 



{n)d 



where <j>{ri) and d are the values of the radial selection function 
and the completeness (0 ^ d ^ 1) for each galaxy, respectively. 
The observed velocity of the Local Group with respect to the CMB 
is given by 



v(r) 



Air 71. 



(8) 



B This is the zero-velocity surface which separates the Local Group from 
the field that is expanding with the Hubble flow (Courteau & Van Den 
Bergh, 1999). 

7 More complicated relations have been suggested for biasing models, ex- 
amples include non-linear and 'stochastic' relations. Also, the halo model 
of clustering involves a different biasing postulation. 
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where n is the mean galaxy density of the survey and f, is the unit 
vector of the galaxy's position. The sum in the equation is over all 
galaxies in the sample that lie in the distance range r m in < < 
Tmax- Calculated this way, the velocity vector does not depend on 
the Hubble constant (h cancels out). 

If the galaxies are assumed to have been drawn by a Poisson 
point process from an underlying density field, then it is straightfor- 
ward to calculate the shot noise errors. The shot noise is estimated 
as the rms of the cumulative variance, a'j n given by 



a ™= E 



(9) 



47rn J t-r 1 \rfcj>(ri)Ci 
The shot noise error per dipole component is am = o" sn /\/3. 

3.2 Flux Weighted Dipole 

For this model, each galaxy is a 'beacon' which represents the 
underlying mass. This is characterised by the mass-to-light ratio 
T = M/L. T is probably not constant and varies with galaxy mor- 
phology (e.g. Lanzoni et al. 2004) but mass-to-light ratios of galax- 
ies vary less in the near-infrared than in the optical (e.g. Cowie et 
al. 1994; Bell & de Jong 2001). In the context of dipole estimation, 
this model of galaxy formation implies that the Newtonian gravita- 
tional acceleration vector for a volume limited sample is 



g( r ) 



= 4ttG 



(10) 



where the sum is over all galaxies in the Universe, (M/L) is the av- 
erage mass-to-light ratio and Si = Li/Aivr 2 is the flux of galaxy i. 
The peculiar velocity vector is derived by substituting Equation llOl 
into the second line of Equation[2] For a flux limited catalogue the 
observed LG velocity is 



v(r) 



87rG/(n m ) 

3£/ofim&L 



77 E W ^ SA 



(ID 



where &l is the luminosity bias factor introduced to account for 
the dark matter haloes not fully represented by 2MRS galaxies and 
WLi is the weight assigned to galaxy i derived in the next section. 
The mass-to-light ratio, assumed as constant, is given by 



/ M 
\ L 



Pm 
PL 



(12) 



87rGpL 

where pl is the luminosity density so Equation llll is rewritten as 



v(r) 



Hp f (Sin 



E 



(13) 



The flux weighting method (originally proposed by Gott) has 
been applied extensively to two-dimensional galaxy catalogues 
(e.g. Yahil, Walker & Rowan-Robinson 1986, Villumsen & Strauss 
1987, Lahav, Rowan-Robinson & Lynden-Bell 1988) and most re- 
cently to 2MASS XSC (Mailer et al. 2003). Since these surveys 
lack radial information, the dipoles were calculated by either as- 
suming = 1 ( e -g- Mailer et al. 2003) or by using a luminos- 
ity function based on a redshift survey in a section of the two- 
dimensional catalogue (e.g. Lahav, Rowan-Robinson & Lynden- 
Bell 1988). In either case, it was not possible to determine the 



convergence of the dipole as a function of redshift. The three- 
dimensional IRAS dipoles (Strauss et al. 1990, Webster, Lahav 
& Fisher 1997, Schmoldt et al. 1999 and Rowan-Robinson et al. 
2000) were derived using the number weighted scheme because the 
IRAS catalogues are biased towards star forming galaxies with 
widely varying mass-to-light ratios resulting in a very broad lu- 
minosity function making it difficult to estimate the distance. The 
2MRS is mainly sensitive to total stellar mass rather than instanta- 
neous star formation rates (e.g. Cole et al. 2001) and consequently 
the 2MRS mass-to-light ratios do not have as much scatter. Thus, 
for the first time, the 2MRS enables the determination the conver- 
gence of the flux dipole as a function of distance. There are many 
advantages to using the flux weighted model for the dipole calcu- 
lation and these will be discussed in the coming sections. 

For the flux weighted case, the weighting function is derived 
as follows: Let pl(L 0) be the luminosity density in the volume 
element SV of a volume limited catalogue. In this case the dipole 
velocity is simply 



v(r) 



E 



(14) 



In practice, however, we have a flux limited catalogue with S ^ 
£ii m so only galaxies with luminosity L ^ L\im = 47rr 2 Sn m are 
included in the survey. Thus the total luminosity in the infinitesimal 
volume SV is 

<5Vp L (L > 0) = iobs + SVp L {L < Uim) (15) 

where L b s = SVp^L ^ Lu m ) is the observed luminosity and 
8Vpl(L < Lii m ) is the luminosity that was not observed due to 
the flux limit of the survey. Substituting 

Lobs 



SV = 

Pl(L ^ Li im ) 

into Equation l 1 5 I vields 



(16) 



SVp^L^O) 



Lob 



1 + 



1 + 



Pl(L < Lli D 
p L (L > Ll in 



p L (L^0)-p L (L^Li in 



p L (L ^ Lu 



Pl(L > 0) 



(17) 



Pl(L ^ Ll im ) V(L ^ Liim) 

where tp(L ^ Ln m ) is the flux weighted selection function. In 
Figure |4] the interpolated fit ip(L Lu m ) for the 2MRS galaxies 
is shown as a function of redshift. 

Thus, the overall weight factor, wl, is 

' (18) 



The luminosity density of the 2MRS is 



= (7.67±1.02)xl0 8 L hMp C - 3 (19) 



where V is the survey volume. The value of pl is in good agree- 
ment with the value derived by Kochanek et al. (2001), pi, = 
(7.14 ± 0.75) x 10 8 L Q /iMpc" 3 . We note that the number 
weighted selection function, </>(r), drops with distance faster than 
the luminosity weighted selection function, lp(r). At large dis- 
tances, we observe only the most luminous galaxies, so the amount 
of 'missing' luminosity from a volume of space is not as big as the 
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Figure 4. The flux weighted selection function as a function of redshift. 

number of 'missing' galaxies. Therefore, as shown below, the flux 
weighted dipole is more robust at large distances than the number 
weighted dipole. 

For the flux weighted scheme, the shot noise is estimated as 

2 / c - \ 2 



E 



(20) 



As the exact shot noise effects for the different models of galaxy 
formation are difficult to model, we will use the Poisson estimate 
above as an indicator of uncertainties (see also Kaiser & Lahav 
1989). However, we note that the quoted uncertainties overestimate 
the noise at small distances where the survey is volume limited and 
underestimate it at large distances where only the brightest galaxies 
are sampled. The uncertainties and the variation in the mass-to-light 
ratios relation also affect the calculation but they are not accounted 
for in this analysis. 

3.3 The Redshift-Space Effects 

It is well known that the peculiar velocities distort the pattern of 
density enhancements in redshift-space. The peculiar acceleration 
of the Local Group calculated using redshifts instead of real dis- 
tances will differ from the actual LG acceleration (Kaiser 1987 and 
Kaiser & Lahav 1989). This effect, referred to as the rocket effect, 
is easily visualised by supposing that only the Local Group has a 
velocity in an homogeneous universe without any peculiar veloci- 
ties. If the LG frame redshifts are used as distance indicators then 
there will be a spurious contribution from the galaxies that are in 
the direction of the LG motion. The prediction for net spurious ac- 
celeration is given by (Kaiser & Lahav 1988): 



fs-pur = -rV(O) 



2+ 



d\n<f) 
d\nr 



'-max, ) ^inax 



(21) 



where z m i n is the minimum, z max is the maximum redshift of the 
survey and z vo i is the redshift for which the survey is volume lim- 
ited (cz vo i ^4500 kms -1 forthe2MRS). 

The rocket effect is very important for the number weighted 
LG dipole calculation because of the dependence on n in Equa- 



tion^ For the 2MRS, the predicted error from the rocket effect is 
a contribution to the total dipole by roughly v apur « — 0.6v(O). 
There are two ways to overcome this error. One is to work in real- 
space instead of redshift-space. This will be discussed in a forth- 
coming paper where the Local Group dipole will be calculated us- 
ing the Wiener reconstructed real-space density field. The other one 
is to use the flux weighted model. For the flux weighted LG dipole, 
the rocket effect is almost negligible as it plays a role only in the 
determination of the radii of the concentric spheres within which 
the dipole is calculated. 



3.4 The Reference Frames 

Brunozzi et al. (1995) and Kocevski, Mullis & Ebeling (2004) 
claim that the dipole in the LG and the CMB frames are over- and 
under-estimates of the real dipole, respectively. The redshift of an 
object is defined by 



cz = Hor + (v(r) - v(0)) ■ f 



(22) 



where v(r) is the peculiar velocity of the object and v(0) is the 
observer's peculiar velocity. In the LG frame |v(0)j = 627 ± 
22kms~ 1 and in the CMB frame |v(0)| = 0kms _1 by definition. 
Therefore, the redshift of a galaxy that has the same direction of 
motion as the LG would be larger in the CMB frame than that in the 
LG frame. In this case, since the acceleration vector is proportional 
to 1/r 2 , the amplitude of the dipole in the LG frame is expected to 
be larger than the amplitude of the dipole in the CMB frame. As the 
dipole is thought to be dominated by the nearby objects that partic- 
ipate together with the LG in a bulk motion (i.e. v(r) ~ v(0) so 
that czlg ~ Hor), it is often assumed that the real LG dipole is 
closer to the dipole in the LG frame than that in the CMB frame. On 
the other hand, Branchini & Plionis (1996) find that the real-space 
reconstruction of the LG dipole gives a result halfway between the 
LG frame and the CMB frame values. 

We perform the analysis using both the LG and the CMB 
frame redshifts. All galaxies are referenced to the rest frame of the 
LG using the transformation in consistency with who use the same 
conversion Courteau & Van Den Bergh (1999): 

czlg = czhei— 79 cos(7) cos (6) +296 sin(7) cos(6)— 36sin(6)(23) 

where z^ei is the heliocentric redshift and I and b are the longitude 
and the latitude of the galaxy in the Galactic coordinates, respec- 
tively. We convert from the LG frame to the CMB frame using 

czcmb = czlg + v LG [sm(b) sm(b LG ) 

+ cos(fe) cos(6lg) cos(\Ilg — l\)], (24) 

where vlg is the amplitude of the LG velocity with respect to the 
CMB and (Ilg, 6lg) is the direction of its motion. We use the 
CMB dipole value of Bennett et al. (2003). Using the first year of 
data from WMAP, they find that the Sun is moving at a speed of 
369.5±3.0 km s~\ towards (I = 263.85° ± 0.10°, b = 48.25° ± 
0.40°). Using the revised values of the motion of the Sun relative 
to the LG a velocity of 306±18 km s _1 towards (I = 99° ± 5°, 
b = -4° ± 4°) derived by Courteau & Van Den Bergh (1999), we 
find a LG velocity relative to the CMB of vlg = 627± 22 km s~ 1 , 
towards (l LG = 273° ± 3°, b LG = 29° ± 3°). 

The choice of reference frames highlights another advantage 
of the flux weighted dipole calculation. As the redshifts only enter 
the calculation in the determination of the radius of the concentric 
spheres, the results are robust to changes in reference frames. 
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4 DIPOLE RESULTS 

The results are presented in Figures l5l7l The top plots of Figures[5j 
|6|show the amplitudes and three spatial components of the accel- 
eration on the Local group (top) and the convergence of the angle 
between the LG dipole and the CMB dipole (bottom) as a function 
of distance in two reference frames. For these plots, the galaxies 
in the masked regions are interpolated from the adjacent regions 
(Method 2). The right panel in each figure shows the results for the 
flux weighted dipole and the left panels show the results for the 
number weighted dipole. Figure|5]is for the the Local Group Frame 
and Figure|6|is in the CMB frame. As discussed in the next section, 
the results for the filling Method 1 where the galaxies are sampled 
randomly do not look very different than the results in Figures [5] 
&[6]and thus are not shown. Figure compares the direction of 
the LG dipole estimate to that of the CMB and other LG dipole 
measurements. 

We give the results for the flux weighted dipole calculated us- 
ing the second method of mask filling in Table|2| Column 1 is the 
radii of the concentric spheres within which the values are calcu- 
lated; columns 2 and 3 are the amplitude of the velocity vector and 
the shot noise divided by Sl m °' 6 /&L, respectively; Columns 4, 5 and 
6 show the direction of the velocity vector and its angle to the CMB 
dipole. The first line gives the results in the LG frame and the sec- 
ond line gives the results in CMB frame. Table[3]is structured in the 
same way as Table |2|however the analysis excludes five galaxies. 

4.1 The Tug of War 

In Figures|5|&|6| the LG velocity is dominated by structure within a 
distance of 60 km s _ 1 (except for the CMB frame number weighted 
dipole where the contribution from the distant structure is over- 
estimated.). The 'tug of war' between the Great Attractor and the 
Perseus-Pisces is clearly evident. The dip in the velocity vector is 
an indication that the local flow towards the Great Attractor 8 is 
counteracted by the Perseus-Pisces complex in the opposite direc- 
tion. If we take out 420 galaxies in the Perseus-Pisces ridge (defined 
by -40 sC b < -10, 110 < b sC 130, 4600 kms -1 sC cz < 6000 
kms ) and recalculate the convergence, the dip almost disappears 
and the convergence is dominated by the Great Attractor. This leads 
us to conclude that the Perseus-Pisces plays a significant role in the 
gravitational acceleration of the LG. 

4.2 Filling the Zone of Avoidance 

The choice of method used to fill the masked regions does not have 
much effect on the results. The convergence of the misalignment 
angle for the second method is slightly more stable than the first 
method and the overall direction of the LG dipole is closer to the 
CMB dipole. Since the Galactic z component is least affected by 
the zone of avoidance, the discrepancy between the amplitudes in 
each plot comes mainly from the Galactic x and y components. 
The direction of the dipole is 2°-3° closer to the CMB vector for 
the second method at distances where the Great Attractor lies. Of 
course, one cannot rule out the possibility that there may be im- 
portant contribution to the dipole from other structures behind the 
zone of avoidance. The Great Attractor is most likely centred on 
the Norma Cluster (Z ^325°, b «-7 Q , Kraan-Korteweg et al. 1996) 

8 By 'Great Attractor', it is meant the entire steradian on the sky centred at 
(I ~ 310°,f> ~ 20°) covering a distance of 20 h' 1 Mpc to 60 h" 1 Mpc. 



which lies very close to the obscured plane. Although, the 2MRS 
samples the Norma cluster much better than the optical surveys, 
the latitude range |6| < 5 which is still obscured in the 2MRS may 
have structure that play an important role in the dipole determina- 
tions. In fact, Scharf et al. (1992) and Lahav et al. (1993) point out 
that there is significant contribution to the local flow by the Pup- 
pis complex at low galactic latitudes that are not sampled by the 
2MRS. 

Kraan-Korteweg & Lahav (2000) point out that since the 
dipole is dominated by local structures, the detection of nearby 
galaxies can be more important to the dipole analyses than the de- 
tection of massive clusters at larger distances. We test this by ex- 
cluding the five most luminous nearby galaxies (Maffei 1 , Maffei 2, 
IC342, Dwingeloo 1 and M81) from our analysis. Remarkably, the 
direction of the resultant dipole moves much closer to that of the 
CMB (see Table|3j. All of these galaxies expect M81 lie very close 
to the Zone of Avoidance and they are excluded from most dipole 
analyses either because they are not in the catalogue (e.g. Rowan- 
Robinson et al. 2000) or they are masked out (e.g. Mailer et al. 
2003). In fact, when we change our mask to match that of Mailer et 
al. (2003) and keep M81 in the analysis our resulting dipole is only 
3° degrees away from the dipole calculated by Mailer et al. (2003). 
This illustrates the importance of the nearby structure behind the 
Zone of Avoidance. 

The comparison of Tables|2|and|3|also highlights the vital role 
non-linear dynamics induced by nearby objects play in dipole cal- 
culations. Mailer et al. (2003) investigate the non-linear contribu- 
tion to the LG dipole by removing the bright galaxies with K s < 8. 
They report that the LG dipole moves to within a few degrees of the 
CMB dipole. They repeat their calculations for the PSCz survey and 
observe the same pattern. We do not observe this behaviour. When 
we remove the objects brighter that K 3 = 8 (428 galaxies), the 
misalignment angle of the resulting dipole decreases by 3° for the 
flux weighted dipole and remains the same for the number weighted 
case. The dipole amplitudes decreases substantially in both cases, 
notably in the case of the flux dipole, suggesting that the brightest 
2MRS galaxies play a significant role in inducing the LG velocity. 



4.3 The Choice of Reference Frames 

The number weighted LG dipole looks very different in differ- 
ent reference frames (Figures [5] and |6j whereas the flux weighted 
dipole is almost unaffected by the change. The number weighted 
dipole is similar to the flux weighted dipole in the LG frame. Thus, 
we conclude that it is more accurate to use the LG frame redshifts 
than that of the CMB frame. 



4.4 The Choice of Weighting Schemes 

The amplitudes of the number weighted LG dipole and the flux 
weighted LG dipole are very similar in the LG frame. However, 
the convergence of the misalignment angles of the flux and the 
number dipoles is quite different in both frames, especially at large 
distances. The angle of the flux weighted dipole is closer to the 
CMB dipole than its number weighted counterpart at all distances. 
With either weighting scheme, the dipoles are closest to the CMB 
dipole at a distance of about 5000 kms -1 and move away from the 
CMB dipole direction further away. However, the change in the di- 
rection of the flux weighted dipole is much smaller compared to 
the number weighted dipole and there is convergence within the 
error bars by 6000 kms -1 . The misalignment angles in the LG 
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Table 1. The values derived for f2 m , j3. 



f2m°' 6 /t | L from the flux weighted scheme 


= 0.40 ±0.09 


f2m°' 6 /b from the number weighted scheme 


= 0.40 ±0.08 



frame at 130 h~ J Mpc are 21° and 37° for the flux and the num- 
ber dipoles, respectively. The discrepancy is probably mainly due 
to the fact the number dipole is plagued with errors due to the lack 
of peculiar velocity information. In fact, when we use just the red- 
shift information instead of the distance measurements (see Section 
2.4) the number dipole moves m 7° towards the CMB dipole. The 
flux dipole assumes that the mass traces light whereas the number 
dipole assumes that all galaxies have the same mass. The former 
assumption is of course more valid, however, since the amplitudes 
of the dipoles are so similar in the LG frame, we conclude that the 
equal mass assumption for the number weighted dipole do not in- 
troduce large errors and that the discrepancy results from the errors 
in distance. 

In all figures, v x and v y change with distance more rapidly in 
the number weighted scheme than the flux weighted. At further dis- 
tances, the flux weighted v x flattens whereas its number weighted 
counterpart continues to grow. It is expected that the (x, y) direc- 
tions are particularly sensitive to the shape of the zone of avoid- 
ance, although it is not obvious why the flux weighted components 
remain so robust. Assuming the dipole has converged we can ob- 
tain values for f2 m °' 6 /6 (number weig hted) and f2 m ° D /6 L (flux 
weighted) by comparing the amplitude of our dipole estimates to 
the CMB dipole. These values are summarised in TableQ The val- 
ues are quoted in the LG frame at 9 13000 kms _1 using the sec- 
ond mask and with the luminosity density value derived earlier, 
p L = 7.67 ± 1.02 x 10 8 L /iMpc" 3 . The errors take the shot 
noise, the uncertainties in the CMB dipole and pl (for the flux lim- 
ited case) into account. The f3 values obtained for the two different 
weighting schemes are in excellent agreement suggesting that the 
dark matter haloes are well sampled by the survey. Our value for /3 
is also in good agreement with results from 2MASS (Pike & Hud- 
son 2005) and IRAS surveys (e.g. Zaroubi et al. 2002, Willick & 
Strauss 1998). In order to calculate the uncertainties introduced by 
the errors in the galaxy redshifts, ten realisations of the 2MRS cata- 
logue are created with each galaxy redshift drawn from a Gaussian 
distribution with standard deviation equal to its error 10 . It is found 
that the scatter in the dipole results due to the errors in redshifts 
are very small compared to the shot noise errors and thus are not 
quoted. 

4.5 The Milky Way Dipole 

We also investigate the acceleration on our galaxy due to other eight 
members of the LG excluded from the LG dipole analysis. As ex- 
pected, the flux dipole is strongly dominated by Andromeda (M31) 
with an amplitude of v/f3 ~ 220kms~ 1 directly towards M31 
(Iw 121.4° , b « —21.7°), confirming that near-infrared fluxes are 
good tracers of mass. The number weighted dipole which assumes 
that the galaxies have the same weight gives a similar amplitude 
of v//3 w 190 km submits direction (I « 104.6°, b « -21.6°) 

9 This is the distance beyond which the shot noise becomes too high (over 
10% for the number weighted analysis). 

10 The mean value of the redshift measurement errors is 30 kms" 1 . 



is skewed towards NGC 6822 (I « 25.3°, b « -18.4°) which 
lies further away from the other seven galaxies that are grouped 
together. 



5 DISCUSSION 

In this paper, we calculate the 2MRS dipole using number and flux 
weighting schemes. The flux weighted dipole bypasses the effects 
of redshift space distortions and the choice of reference frames giv- 
ing very robust results. 

Our dipole estimates are dominated by the tug of war be- 
tween the Great Attractor and the Perseus-Pisces superclusters 
and seemingly converge by 6000 kms -1 . The contribution from 
structure beyond these distances is negligible. The direction of 
the flux dipole (l=251°±12°,b=37°± 10 °) is in good agree- 
ment with the 2MASS dipole derived by Mailer et al (2003) 
(1=264.5° ±2° ,b=43.5°± 4 °). The difference in results is prob- 
ably due to the fact that they use a higher latitude cutoff in the 
mask (|6| < 7°) and exclude all galaxies below this latitude. We 
confirm this by changing our treatment of the Zone of Avoidance 
to match theirs. We find that the flux dipole is very close to their 
dipole direction. Their limiting Kron magnitude is K s — 13.57 
which corresponds to an effective depth of 200 h^ 1 Mpc. As their 
sample is deep enough to pick out galaxies in the Shapley Super- 
cluster, the comparison of their dipole value with our values sug- 
gests that the contribution to the LG dipole from structure further 
away than the maximum distance of our analysis is not signifi- 
cant. Following Mailer et al. (2003), when we adopt £3 m = 0.27 
we get 6_l = 1.14 ± 0.25, in good agreement with their value of 
b — 1.06 ± 0.17 suggesting that the 2MRS galaxies are unbiased. 
We note that the 2MRS value for the linear bias is somewhat lower 
than expected considering that the 2MRS has a high proportion of 
early type galaxies which are known to reside mostly in high den- 
sity regions (e.g. Norberg et al. 2001, Zehavi et al. 2002). The val- 
ues we derive for /3 and fl m are consistent with the concordance 
A-CDM model values given their error bars. 

Figure 3 shows that the 2MRS samples the Great Attractor 
region better than the PSCz survey but that the PSCz redshift dis- 
tribution has a longer redshift tail than the 2MRS. Nevertheless, 
PSCz dipole agrees well with that of the 2MRS. Rowan-Robinson 
et al. (2000) derive a value for (3 = 0.75t o ; o £ which is higher 
than the value we derive in this paper. The PSCz sample is biased 
towards star-forming galaxies and thus under-samples the ellipti- 
cals which lay in high density regions. Contrarily, the 2MRS is bi- 
ased towards early-type galaxies. This difference may be the reason 
why they get a higher value for f3. The flux weighted dipole is in 
excellent agreement with the IRAS 1.2 Jy dipole (Webster, La- 
hav & Fisher 1997) which was obtained using a number weighted 
scheme but with an added filter that mitigates the shot noise and de- 
convolves redshift distortions. Our number weighted dipole differs 
from their results. This is probably due to fact that the 2MRS num- 
ber weighted dipole is plagued wit redshift distortions. Webster, 
Lahav & Fisher (1997) obtain the real-space density field from that 
in the redshift-space using a Wiener Filter. In a forthcoming paper, 
we will use the same technique to address this issue. Similarly, the 
2MRS number weighted dipole also differs from the PSCz dipole 
which was calculated using a flow model for massive clusters. 

The analysis of the IRAS QDOT galaxies combined with 
Abell clusters (Plionis, Coles & Catelan 1993) and the X-Ray clus- 
ter only dipole (Kocevski, Mullis & Ebeling 2004 and Kocevski 

6 Ebeling 2005) imply a significant contribution to the LG veloc- 
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Table 2. The convergence values calculated for flux weighted dipole using the second method for the masked region. The columns give (i)Radius of the 
concentric spheres within which the values are calculated; (ii) and (iii) amplitude of the velocity vector and the shot noise divided by /(f2 m ), respectively; 
(iv), (v) and (vi) direction of the velocity vector and its angle to the CMB dipole. The results in the first and second line are given in LG and CMB frames, 
respectively. 



Dist 

i — 1 
km s 


v tot b L /f(nm) 

1 —1 
km s 1 


1 

deg 


b 

deg 


se 

deg 


1000 


590 ±294 
362 ±289 


259 ± 59 
331°±171 ° 


39 ±20 

23° ±23 ° 


42 ±20 

99° ±26° 


2000 


1260 ±318 
993 ±316 


249 ± 17 
222° ±22° 


j i O 1 1 O O 

41 ±13 
41°±15 ° 


25 ±11 
44° ±15 ° 


3000 


1633 ±322 
1334 ±320 


260 ± 12 
241°± 17 ° 


/IaO 1 1 f\ O 

40 ±10 

44°±13 ° 


17 ± 8 
31°±11 ° 


4000 


1784 ±323 
1513 ±322 


264 ± 10 

252° ± 14 ° 


1 AO 1 f\ O 

39 ± 9 
41°±11 ° 


14° ± 7 ° 
22°±9 ° 


5000 


1838 ±324 
1497 ±323 


265 ± 10 

252° ± 13 ° 


36 ± 9 
39°±11 ° 


12 ± 7 

22°±9° 


6000 


1 /TO n 1 O'lC 

1633 ±325 
1438 ±324 


259 ± lU 
250° ± 12 ° 


n a O i r» O 

34 ± 9 
36°±11 ° 


15°± 8 ° 
22°±9° 


7000 


1682 ±325 
1503 ±324 


r\ C f O 1 1 A O 

256 ± 10 

247° ± 12 ° 


36 ± 9 
36°±10° 


18 ± 8 

24°±9° 


8000 


1691 ±326 
1566 ±325 


255 ± ll 

248° ± 12 ° 


TnO l 1 *~\ O 

37 ±10 
39° ±10° 


18 ± 8 
24°±9° 


9000 


1683 ±326 
1573 ±325 


r\ c C O 1 11 O 

255 ± ll 

248° ± 12 ° 


o o O l 1 f\ O 

38 ±10 
39° ±10° 


I aH i o O 

19 ± 8 
24°±9° 


10000 


1674 ±326 
1599 ±325 


^ ^ 1 O 1 11 O 

253 ± ll 
246° ± 12 ° 


o o O l 1 /A O 

38 ±10 
39°±10° 


r\AO 1 o O 

20 ± 8 
26°± 8 ° 


11000 


1677 ±326 
1624 ±325 


O 1 11 O 

253 ± ll 
246° ± 12 ° 


o o O l 1 /"\ O 

38 ±10 
40°±10° 


^ 1 O I o O 

21 ± 8 

26°± 8 ° 


12000 


1676 ±326 
1624 ±325 


c n O 1 11 O 

253 ± ll 
246° ± 12 ° 


^> o O 1 1 /A O 

38 ±10 
40°±10° 


1 O I Ci O 

21 ± 8 
26°± 8 ° 


13000 


1652 ±326 
1629 ±325 


^ r 1 O I 11 O 

25 1 ± ll 

245°± 12 ° 


38 ±10 
39°±10° 


^> 1 O I o O 

21 ± 8 

26°±8° 


14000 


1659 ±327 
1636 ±326 


c 1 O I 11 O 

25 1 ± ll 
245°± ll ° 


TOO I i /\ O 

38 ±10 
39°±10° 


22 ± 8 
26°± 8 ° 


15000 


1640 ±327 
1633 ±326 


c 1 O I 11 O 

25 1 ± ll 
246°± ll ° 


TnO l 1 f\ O 

37 ±10 
39°±10° 


1 O I Ci O 

21 ± 8 

25°± 8 ° 


16000 


163s ±327 
1643 ±326 


251 ±11 
247°± 11 ° 


37 ±10 
38°±10° 


Tl O _|_ O 

21 ± 8 
25°± 8 ° 


17000 


1630 ±327 
1643 ±326 


251°± 11 ° 
247°± 11 ° 


37°±10° 
38°±10° 


21°± 8 ° 
25°± 8 ° 


18000 


1604 ±328 
1631 ±327 


251°± 11 ° 
247°± 11 ° 


37°±10° 
38°±10° 


21°± 8 ° 
25°± 8 ° 


19000 


1591 ±328 
1629 ±327 


251°± 11 ° 
247°± 11 ° 


37°±10° 
38°±10° 


21°± 8 ° 
24°± 8 ° 


20000 


1577 ±328 
1620 ±327 


251°± 12° 
247°± 11 ° 


37°±10° 
37°±10° 


21°± 8 ° 
24°± 8 ° 
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Table 3. Same as Tablel2*lbut the analysis excludes the galaxies: Maffei 1, Maffei 2, Dwingeloo 1, IC342 and M81. 



Dist 
km s _1 


VtotDL / f (Urn) 
km s~ 1 


1 

deg 


b 
deg 


Ou 
deg 


1000 


585 ±254 
181 ±228 


280° ± 30 ° 
360° ±108 ° 


34° ±17 ° 
24° ±29 ° 


25°±9° 
67° ±37 ° 


2000 


1258 ±282 
963 ±261 


263°± 13 ° 
259°± 18 ° 


38°±11 ° 
39°±13 ° 


15°±9° 
19° ±12 ° 


3000 


1661 ±285 
1352 ±266 


269° ±9 ° 
265° ± 12 ° 


37°±9° 
41°±10° 


11°±5 ° 
16°±8° 


4000 


1824 ±287 
1571 ±268 


272° ± 8 ° 
269° ±9 ° 


36°±8° 
37°±9 ° 


10°±4° 
12°±5 ° 


5000 


1891 ±288 
1559 ±270 


272° ± 8 ° 
268°±9° 


33°±7 ° 
35°±8° 


8°± 3 ° 
10°±5 ° 


6000 


1678 ±289 
1504 ±271 


268°±8° 
266° ±9 ° 


31°± 8 ° 
32°± 8 ° 


9°±4° 
10°±5 ° 


7000 


1713 ±289 
1551 ±272 


264° ± 8 ° 
263°±9° 


33°± 8 ° 
33°± 8 ° 


11°±6° 
12°±6° 


8000 


1721 ±290 
1611 ±272 


264°±9° 
264°±9° 


35°±8° 
36°± 8 ° 


11°±7° 
12°±7 ° 


9000 


1704 ±290 
1614 ±272 


264°±9° 
264°±9° 


35°± 8 ° 
36°± 8 ° 


12°±7 ° 
12°±7 ° 


10000 


1691 ±290 
1634 ±272 


262°±9° 
262°±9° 


36°± 8 ° 
36°± 8 ° 


13°±7 ° 
13°±7 ° 


11000 


1691 ±290 
1658 ±273 


262°±9° 
262°±9° 


36°± 8 ° 
37°± 8 ° 


13°±7 ° 
14°±7 ° 


12000 


1690 ±291 
1658 ±273 


262°±9° 
261°±9° 


35°± 8 ° 
37°± 8 ° 


13°±7 ° 
14°±7 ° 


13000 


1665 ±291 
1663 ±273 


261°±9° 
261°±9° 


35°±9 ° 
36°± 8 ° 


14°±7 ° 
14°±7 ° 


14000 


1671 ±291 
1668 ±273 


260°±9° 
260° ±9 ° 


35°± 8 ° 
36°± 8 ° 


14°±7 ° 
14°±6° 


15000 


1652 ±291 
1670 ±274 


260° ±9 ° 
261°±9° 


35°±9 ° 
36°±8° 


14°±7 ° 
13°±6° 


16000 


1653 ±292 
1683 ±274 


261°±9° 
261°± 8 ° 


34°±9° 
35°± 8 ° 


13°±7 ° 
13°±6° 


17000 


1647 ±292 
1684 ±274 


261°±9° 
261°± 8 ° 


34°± 8 ° 
35°±8° 


13°±7 ° 
13°±6° 


18000 


1619 ±292 
1672 ±274 


260°±9° 
261°±9° 


34°±9° 
35°± 8 ° 


14°±7 ° 
13°±6° 


19000 


1606 ±292 
1672 ±274 


260°±9° 
262° ± 8 ° 


34°±9° 
35°± 8 ° 


14°±7 ° 
13°±6° 


20000 


1594 ±293 
1665 ±275 


261°±9° 
262° ± 8 ° 


34°±9° 
34°± 8 ° 


13°±7 ° 
12°±6° 
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Figure 5. Top: Three components and the magnitudes of the acceleration of the Local Group due to galaxies within a series of successively larger concentric 
spheres centred on the local group in the Local Group frame. The galaxies in the masked regions are interpolated from the adjacent regions (Method 2). 
Left panel is the number weighted velocity and right panel is the flux weighted velocity. The growth of the estimated shot noise is also shown. Bottom: 
Convergence of the direction of the LG dipole where the misalignment angle is between the LG and the CMB dipoles. The dotted lines denote Icr errors from 
shot noise. Left plot is the direction of the number weighted LG dipole and right plot is the direction of the flux weighted LG dipole. We note for the number 
weighted dipole the dramatic increase in shot noise beyond 15000 kms -1 , where the dipole's behaviour cannot be interpreted reliably. 
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Figure 6. Same as in Figure|5]but in CMB frame. 
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Figure 7. The triangles show the dipole directions at 1 3000 km s 1 derived in this paper. The direction of the number weighted 2MRS dipole in the LG frame 
is shown in magenta (LG-N.W., I = 231°, b = 42°). The red triangle shows the direction of the flux weighted 2MRS in the LG frame (LG-F.W., I = 251°, 
b = 38°). The blue triangle shows the direction of the flux weighted 2MRS in the LG frame excluding Maffei 1, Maffei 2, M81, IC342 and Dwingeloo 1 
(LG-F.W., I = 261°, b = 35°). The red star shows the CMB dipole direction (2 = 273°, b = 29°). The green triangle is the number weighted LG dipole in 
the CMB frame (CMB-N.W., I = 218°, b = 33°); the yellow triangle is the flux weighted LG dipole in CMB frame (CMB-FW., I = 245°, b = 39°). The 
green circle is the 2MASS dipole (/ = 264°, b = 43°, Mailer et al. 2003). The magenta square is the IRAS 1.2-Jy dipole (I = 247°, b = 37°, Webster 
et al. 1997). The blue upside-down triangle is the IRAS PSCz dipole (I = 253°, b = 26°, Rowan-Robinson et al. 2000). Contours are drawn at constant 
misalignment angles. Also shown are the Virgo Cluster (I = 280°, b = 75°), the Hydra Cluster (I = 270°, b = 27°), the Centaurus Cluster (/ = 302°, 
b = 22°) and A3558 (/ = 312°, b = 31°). 



ity by the Shapley Supercluster. Kocevski & Ebeling (2005) report 
a significant contribution to the LG dipole (56%) from distances 
beyond 60 h^ 1 Mpc. The discrepancy between their results and 
ours is possibly due to the fact that the 2MRS is a better tracer of 
the galaxies at nearby distances, whereas the X-ray cluster data are 
better samplers of the matter distribution beyond 150 h^ 1 Mpc. 

The misalignment angle between the LG and the CMB dipole 
is smallest at 5000 km s _1 where it drops to 12° ±7° and increases 
slightly at larger distances presumably due to shot-noise. This be- 
haviour is also observed in the other dipole analyses (e.g. Web- 
ster, Lahav, Fisher 1997 & Rowan-Robinson et al. 2000). This is a 
strong indication that most of the LG velocity is due to the Great 
Attractor and the Perseus-Pisces superclusters. Of course, we still 
cannot rule out a significant contribution from Shapley as we do 
not sample that far. However, it may be more important ask what 
the velocity field in the Great Attractor region is. In other words, 
whether we observe a significant backside infall towards the Great 
Attractor. 



The smallest misalignment angle at 13000 kms _1 is 21°±8°, 
found for the LG frame, flux weighted scheme using the second 
mask. This misalignment can be due to several effects: 

• The analysis uses linear perturbation theory which is correct 
only to first order 8. There may be contributions to the LG dipole 
from small scales which would cause gravity and the velocity vec- 
tors to misalign, even with perfect sampling. However, Ciecielag, 
Chodorowski & Kudlicki (2001) show that these non-linear effects 
cause only small misalignments for the IRAS PSCz survey. How- 
ever, removing the five most luminous nearby galaxies moves the 
flux dipole to (I = 261° ±9°, fa = 34° ±9°,cz=20000 kms" 1 ), 8° 
closer to that of the CMB. This suggests that the non-linear effects 
might be very important in dipole determinations. 

• The sampling is not perfect and the selection effects of the 
surve will increase the shot noise-errors especially at large dis- 
tances causing misalignments. 

• There may be uncertainties in the assumptions in galaxy for- 
mation and clustering. For example the mass-to-light ratios might 
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differ according to type and/or vary with luminosity or the galaxy 
biasing might be non-linear and/or scale dependent. 

• There may be a significant contribution to the LG dipole from 
structure further away than the maximum distance of our analysis. 

• The direction of the LG dipole may be affected by nearby 
galaxies at low latitudes which are not sampled by the 2MRS. In 
the future, the masked regions will be filled by galaxies from other 
surveys such as the ongoing HI Parkes Deep Zone of Avoidance 
Survey (Henning et al. 2004) as well as the galaxies that are sam- 
pled from the 2MRS itself. 

Our initial calculations of the expected LG acceleration (c.f. Lahav, 
Kaiser & Hoffman 1990; Juszkiewicz, Vittorio & Wyse 1990) sug- 
gest that the misalignment of 21°±8° is within 1 a of the dipole 
probability distribution in a CDM Universe with Q m = 0.3. In 
a forthcoming paper, the cosmological parameters will be con- 
strained more vigourously using a maximum likelihood analysis 
based on the spherical harmonics expansion (e.g. Fisher, Scharf & 
Lahav 1994; Heavens & Taylor 1995) of the 2MRS density field. 
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